Normal stresses at the gelation transition 
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A simple Rouse-type model, generalised to incorporate the effects of chemical crosslinks, is used 
to obtain a theoretical prediction for the critical behaviour of the normal-stress coefficients \Pi and 
\p2 at the gelation transition. While the exact calculation shows \&2 = 0, a typical result for these 
types of models, an additional scaling ansatz is used to demonstrate that ^i diverges with a critical 
exponent t — k + z. Here, k denotes the critical exponent of the shear viscosity and z the exponent 
governing the divergence of the time scale in the Kohlrausch decay of the shear-stress relaxation 
function. For crosslinks distributed according to mean-field percolation, this scaling relation yields 
£ — 3, in a accordance with an exact expression for the first normal-stress coefficient based on a 
replica calculation. Alternatively, using three-dimensional percolation for the crosslink ensemble we 
find the value £ « 4.9. Results on time-dependent normal-stress response are also presented. 

PACS numbers: 64.60.Ht, 61.25.Hq, 61.20.Lc 



I. INTRODUCTION 

Chemical gelation is the process of randomly introduc- 
ing crosslinks between the constitutents of a (macro-) 
molecular fluid. One way to investigate the effects of the 
crosslinks on the fluid dynamics consists in measuring the 
stresses the crosslinked fluid builds up when subjected to 
a simple shear flow. For an incompressible, isotropic fluid 
one can experimentally access m three independent com- 
ponents of the stress tensor a: the shear stress a xy and 
the first and second normal stress differences a xx — a yy 
and a yy — a zz . For static shear flows these give rise to 
three independent material functions: the shear viscos- 
ity 77 and the first and second normal-stress coefficients 
"Fi and ^2- Generally speaking, both Newtonian and 
non-Newtonian fluids possess a non-vanishing shear vis- 
cosity. But, whereas for a Newtonian fluid both 'Fi and 
^2 are always zero, it is precisely the non-vanishing of 
'Fi that explains a number of characteristic effects known 
for e.g. polymeric liquids g], see also §2.3 in [[[J. On the 
other hand, even for non-Newtonian fluids ^2 is typi- 
cally found to be very small as compared to \&i, and the 
"Weissenberg hypothesis" , ^2 = 0, is a good approxima- 
tion in these cases ||. It also seems that ^2 is not as 
well investigated experimentally as \&i. 

In the context of gelation one is particularly interested 
in the dependence of these stresses on the crosslink con- 
centration c. Universal critical behaviour is expected 
to occur at the gelation transition, that is, at the crit- 
ical concentration c cr ;t, where the fluid (sol) undergoes 
a sudden phase change into an amorphous solid state 
(gel). As far as shear stress is concerned, there exist nu- 
merous experimental investigations on the static shear 
viscosity and on the time-dependent shear-stress relax- 
ation function. The experimentally measured values for 
the critical exponent k, which governs the algebraic di- 



vergence of the shear viscosity when approaching c cr it 
from the sol side, scatter considerably and are found in 
the range fc « 0.6 . . . 1.7, see e.g. § § | | || §. The 
origin of this wide spreading is controversially discussed 
and eventually unclear. From a theoretical point of view 
there exists a bunch of competing and partially contra- 
dicting scaling relations which express k in terms of per- 
colation exponents. Each of them relies on heuristic ar- 
guments whose validity is mostly unclear. We refer the 
reader to P, pll for a summary and references. Here we 
only mention the scaling relation k = 2v — f3 which was 
first proposed by de Gennes [fill and rederived by many 
others. Erroneously, it is generally referred to as the 
"Rouse expression" for the viscosity exponent. Here, v is 
the exponent governing the divergence of the correlation 
length and j3 is associated to the gel fraction. For three- 
dimensional bond percolation one would get the value 
(2v — /3)| d=3 ~ 1-35. Recently, the viscosity was exactly 
determined within the Rouse model for gelation in fi [Kj . 
The analysis disproves the above result and shows that 



k = 4>- (3 



(1) 



"Deceased 



is the true scaling relation valid for Rouse dynamics. 
Here, <f> denotes the first crossover exponent of a cor- 
responding random resistor network |l2j, |l3J. When 
inserting fi HG] high-precision data for 4> an d P ob- 
tained from three-dimensional bond percolation, the true 
Rouse value of the viscosity exponent turns out to be 
(<t> — /3) I d=3 ~ 0.71 and agrees with simulations Q on 
a similar model. The discrepancy to de Gennes' result 
above can be attributed to the neglect of the multiir&ct&l 
nature of percolation clusters in [|ll| . Amazingly, the true 
Rouse value k\, w 0.71 differs only little from that 
of another proposal, k — s, by de Gennes |15f| , where 
he alluded to an analogy to the conductivity exponent 
s sa 0.73 of an electrical network consisting of a ran- 
dom mixture of superconductors and normal conductors. 
This close agreement, however, seems to be coincidental. 
In contrast, the authors are not aware of any experi- 
mental or theoretical studies concerning the dependence 



on the crosslink concentration c of normal stresses near 
the gelation transition. This seems all the more surpris- 
ing since there exist many experiments (l^, [l7], [l8[ [ll| on 
both the shear-rate dependence of normal stresses in en- 
tangled or (temporarily) crosslinked polymeric liquids in 
order to explain shear-thinning or shear-thickening phe- 
nomena and on the time dependence of the normal-stress 
response to particular shapes of shear strain. Theoret- 
ical explanations of these experimental findings mainly 
rely on the analysis of transient network models, see e.g. 

Even though Rouse-type models incorporate no other 
physical interactions between monomers apart from con- 
nectivity, they serve as a standard theoretical reference 
in terms of which experimental data are frequently inter- 
preted. Therefore it is important to test their predictions 
as accurately as possible. It is the purpose of this Pa- 
per to use the same generalised Rouse-type model as in 
0, [1(1 p4| to predict the critical behaviour of the normal- 
stress coefficients v&i and ^ at the gelation transition. 
Within this model it will turn out that $2 vanishes for 
all c and that ^i diverges with a critical exponent 



I = k + z 



(2) 



when approaching c cr j t from the sol side. Here, z de- 
notes the exponent governing the divergence of the time 
scale in the Kohlrausch decay of the shear-stress relax- 
ation function. For crosslinks distributed according to 
mean- field percolation (also called "classical theory"), 
this scaling relation yields I = 3, in a accordance with 
an exact expression for ^i based on a replica calcula- 
tion. Alternatively (and more realistically), using three- 
dimensional percolation for the crosslink ensemble we 
find the value i « 4.9. Thus, the model predicts a 
much more pronounced divergence of ^i as compared 
to r\ so that \I/i may serve as a sensitive indicator for 
the gelation transition. We also derive results on the 
time-dependent normal-stress response. In particular, 
the Lodge-Meissner rule, see e.g. §3.4.e in Q, is shown to 
hold for normal-stress relaxation after a sudden shearing 
displacement. 

We hope that these theoretical investigations moti- 
vate corresponding experimental work in order to develop 
more insight on normal stresses in (near) critical gels. 



II. MODEL 

We follow a semi-microscopic approach to gelation 
based on a Rouse-type model for N monomers. The 
monomers are treated as point particles with positions 
Rj(t), i = 1, . . . , N, in three-dimensional space. The mo- 
tion of the monomers is constrained by M randomly cho- 
sen, harmonic crosslinks which connect the pairs (i e ,i' e ), 
e = 1, . . . , M, of monomers and give rise to the potential 
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FIG. 1: Homogeneous linear shear flow (g) 
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Here, the fixed length a > models the overall inverse 
coupling strength, whereas the individual coupling con- 
stants A e are chosen at random. Quite often, only the 
special case A e = 1 has been considered previously. The 
second equality in (Q) introduces the random N x TV- 
connectivity matrix, which encodes all properties of a 
given crosslink realisation. 

Following 1 25, ^6| £7] we employ a simple relaxational 
dynamics 

BIT 

C[ftfl?(*)-«2 rt (Ri(t),*)]=-^s(*) + ^(t) (4) 

without inertial term to describe the motion of the mo- 
nomers in the externally applied velocity field 



t (r,i) := 5 ax K(t)y 



(5) 



with a time-dependent shear rate K(t), see also Fig. [l| 
Here, Greek indices label Cartesian coordinates x,y or 
z. A friction force with friction constant £ applies if the 
velocity of a monomer deviates from the externally ap- 
plied flow field. The crosslinks exert a force —dU/dRi 
on the monomers, in addition to a random, fluctuating 
thermal-noise force obeying Gaussian statistics with zero 
mean and covariance (£"(i)£f (£')) = 2C$ap{>ij5(t — £')• 
Note that we have chosen units in which the inverse tem- 
perature is equal to one. Given the shear flow (||), the 
equation of motion (0) is linear and can be solved exactly 
for each realisation of the thermal noise |1(J . 

To complete the description of the dynamic model, we 
have to specify the statistical ensemble which determines 
the realisations of the crosslinks. We will distinguish two 
cases: 

(i) Mean-field percolation (also called "classical the- 
ory"): each pair of monomers is chosen indepen- 
dently with equal probability M/N 2 , irrespectively 
of the monomer positions in space. As a function 
of the crosslink concentration c := M/N , the sys- 
tem undergoes a percolation transition at a criti- 
cal concentration c cr it = h- For c < c CT a there is 



no macroscopic cluster, and almost all clusters are 
trees @ ||. 

(ii) Three-dimensional bond percolation |3(J, |3l] . 

For either case we assume the random coupling constants 
A e to be distributed independently of the crosslink config- 
uration, as well as independently of each other with the 
same (smooth) probability distribution p(X). Moreover, 
sufficiently high inverse moments 



from arbitrary macromolecules, as is the case in vulcan- 
ization, for example. This general universality argument 
was confirmed |H| by explicit computations of the criti- 
cal behaviour of the shear viscosity with the mean-field 
distribution of crosslinks. 



III. STRESS TENSOR AND NORMAL-STRESS 
COEFFICIENTS 



P« 



dAA~>(A) 



(0) 



of A e shall exist. 

The combined average over crosslink configurations 
and random coupling constants will be denoted by an 
overbar • . Using this notation, we implicitly assume 
that the macroscopic limit N — > oo, M — > oo, M/N — ► c 
is carried out, too. 

Before turning to the analysis of the model we would 
like to comment on the fact that it describes the ran- 
dom crosslinking of single monomers rather than of pre- 
formed polymers. However, this does not mean that the 
applicability of the model is limited to the description 
of random networks built up by polycondensation from 
small structural units. Indeed, we expect from universal- 
ity that details at small length scales are irrelevant for 
the true critical behaviour at the gelation transition so 
that these results will also hold for random networks built 



Due to the externally applied shear flow the crosslinks 
exert shear stress on the polymer system, whose tensor 
components are given in terms of a force-position corre- 
lation lEl |2l 



Oap{t) = 



lim — 



N 



to- 



" =1 



dU 



OR 



3 (*)flf(t) 



(7) 



Here, po denotes the density of monomers. In (|7|) one has 
to insert the explicitly known fin] solutions R^ (£) of the 
Rouse equation (Q) at time t with initial data Ri(io) a t 
time to. In order to ensure that the thermal- noise average 
allows for the description of a possible stationary state of 
the system at finite times t, the time evolution is chosen 
to start in the infinite past, to — > —00, thereby losing all 
transient effects which stem from the initial data. This 
yields m(x for the stress tensor 



J 



<r(t) =X(0)1 



dt'x(t -*')«(*') 



2/ t ,ds«;(s) 1 

" 1 





(8) 



where 1 denotes the 3 x 3-unit matrix and the stress 
relaxation function is given by 

6i 

Co 5 



X (t) := £ Tr 

AV J N 



(1-Eo) cxp(-— r 



(9) 



The symbol Tr in (||) stands for the trace over N x TV- 
matrices, and Eq denotes the projector on the space of 
zero eigenvalues of T, which correspond to translations of 
whole clusters. The associated eigenvectors are constant 
within each cluster and zero outside |l(J, Q. Within 
the simple Rouse model the zero eigenvalues do not con- 
tribute to shear relaxation because there is no force act- 
ing between different clusters. The only contribution to 
stress relaxation is due to deformations of the clusters. 

For a time-independent shear rate n(t) = k it is cus- 
tomary to define a first and second normal-stress coeffi- 
cient by 



*i 



*9 



POK 



&vv O ' z 



Pqk 



(10) 



One deduces immediately from (||) that 

*2=0, 



(11) 



a characteristic result for Rouse-type models. In con- 
trast, the first normal-stress coefficient ^i is non-zero 



*i 



2 V 3 J 



2 1 

N 



Tr 



1-Ep 

r 2 



(12) 



and independent of the shear rate k. 

For a macroscopic system VPi is expected to be a self- 
averaging quantity. Therefore we will calculate the dis- 
order average • of ( |l2| ) over all crosslink realisations 
and all crosslink strengths. To do so it is convenient to 
introduce the averaged density 



D(ry) 



1tv[(1-E )6^-T)] 



(13) 



of non-zero eigenvalues of T. Physically, D describes the 
distribution of relaxation rates in the network (in units 



of £a 2 /6), as is evident from the representation 



and hence 



X(t) = A) 



| °d 7 exp{-^}^( 7 ) (14) 



of the disorder average of the stress relaxation function 
(||). Various properties of the eigenvalue density D are 
discussed in detail in |p4[ . The average \Pi now appears 
as the second inverse moment of D, 



*i 



\mi?-^ 



T 



(15) 



while the disorder- averaged static shear viscosity rj := 
&xy/{poK) is determined fi FlG] by the first inverse mo- 
ment 

Po Jo 6 J 7 

At this point one can already see that \&i serves as 
a sensitive indicator for the gelation transition. Indeed, 
the Jensen inequality |33l implies 



#i > 



2rJ : 



Jo°° d 7£(7) 



>2rj 2 , 



(17) 



> 2fc 



(18) 



with t, respectively k, denoting the critical exponent of 
*i ~ (ccrit - c)~ l , respectively r] ~ (c cr it - c)~ fe - 

In the two following sections we will determine the 
precise Rouse value of £ for the two different types of 
crosslink ensembles described above. 



IV. MEAN-FIELD PERCOLATION 



For mean-field random graphs (i) the second inverse 
moment of the eigenvalue density D was calculated in 
Eq. (38) of |24|] with the help of a replica approach. This 
gives rise to the exact result 



_ 1 ( C,a 2 

*i = - - — 
2 V 3 



8c d - 6c 2 - 5c + 1 
30c(l - 2c) 3 



Pf 



Ac 2 - 3c - 1 



24c(l 
I 



2c) 2 



Pa 



5P 2 - AP 2 
240c 2 



ln(l - 2c) 



(19) 



which is valid for all < c < c cr i t = h. The inverse 
moments P n were defined in (p). From (113) we read off 
the critical divergence 

_ /fn 2 \? p 2 

*i~(y) ^e~\ e:=c clit -cl0 (20) 

at the gelation transition, and hence the critical exponent 



3. 



For c — ► one expands 



*i 



(l)i-*)' 



(21) 



(22) 



Fig. H displays <3/i in units of (£a 2 /3) 2 as a function of c 
for the special case Pi = P2 = 1. 

It is the merit of the mean-field percolation ensemble 
that it allows for a variety of exact analytical calculations. 
However, since the probability for a crosslink to occur 
does not depend on the monomers' positions in space, 
this ensemble is believed to provide a fairly unrealistic 
description for three-dimensional gels. For this reason 
we consider an alternative crosslink ensemble in the next 



section, which has been successfully used 34 to explain 
static properties of polymer systems. 



V. THREE-DIMENSIONAL BOND 
PERCOLATION 

For this ensemble of crosslinks the second inverse mo- 
note that in this 



ment (15) of the eigenvalue density D, 
section we emphasize the dependence on e :— c cr it — c in 
the notation of various quantities - is not known analyt- 
ically. In order to proceed we assume that D e follows a 
scaling law 



De{l)~l A - l f(l*{e)h) 



(23) 



It 



close to the critical point and for small enough 7 
is determined by a typical relaxation rate 7*(e) ~ e z , 
which vanishes when approaching the critical point, and 
a scaling function fix) that tends to a non-zero constant 
for x — > and decays faster than any inverse polynomial 
for x — > 00. In particular, this gives the power-law be- 
haviour D E= o(-f) ~ 7 A_1 asymptotically for 7 — > at 
criticality, in agreement with experiments |35| , p6| . The 
measured exponent values, however, scatter considerably, 
A w 0.4 . . . 0.8, and seem to depend on the mass of the 
crosslinked molecules p?J . Note that on general grounds 
the exponent A has to be positive, because otherwise 
^£=0(7) would not be integrable at 7 = 0, in contradic- 




0.125 0.25 0.375 0.5 



FIG. 2: First normal-stress coefficient ( |L9[ ) in units of (£a 2 /3) 2 
as a function of c for P 1 = P 2 = 1. 



tion to the definition (|13|). The scaling law ( j23| ) yields, 
via the Laplace transform (O), the scaling law 



when using three-dimensional bond percolation to gen- 
erate the crosslink ensemble. Concerning A, we are only 
aware of [M , where this exponent is determined for the 
Rouse model at hand without any further assumptions. 
It was done by numerical computations of the eigenvalue 
density ( |l3|) and yields A rs 0.83. But, as compared to 
k, we suspect the numerical accuracy of the result to be 
rather poor. Yet, using these values, Eq. (|27|) predicts 



£«4.9 



(28) 



Xe(t)-e zA g(t/r*(s)) 



(24) 



for the exponent of the first normal-stress coefficient \&i. 
If, instead, one ignored the mw^zfractal structure of per- 
colation clusters in employing the wrong scaling relations 
k = 1v — (3 and t — dis, where t is the critical exponent 
of the elastic modulus in the gel phase, one would arrive 
at the value A « 0.66. This would yield the considerably 
lower result t « 2.8. Thus, it is of importance to im- 
prove the accuracy of the exact numerical computation 
of A within the Rouse model. 

Finally, we would like to point out that for mean-field 
percolation the scaling relation (£7|) is consistent with the 
exact result presented in the preceding section. For, in 
this case the model yields k — |j, 10 1, z — 3 and A = 1 
p4f, and thus (p7|) gives I — 3 in accordance with (|2l|). 



for the long-time behaviour of the stress relaxation func- 
tion. Here, the scaling function obeys g(x) ~ x~ A 
for x — > and the typical relaxation time T*(e) := 
Ca 2 /[67*(e)] ~ e~ z diverges when approaching the crit- 
ical point. Precisely at the critical point one finds an 
algebraic long-time decay X e= o(4) ~ t~ A . Dynamical 
scaling relates A to z and to the exponent k of the shear 
viscosity 



A = (z-k)/z, 



(25) 



see e.g. [£4|, [37|. For x — + oo the scaling function g(x) has 
to decay like a stretched exponential in order to accom- 
modate the experimentally found pll Kohlrausch decay 



X e>0 (i) ~ exp{-[t/r*(e)]«} 



(26) 



of the stress-relaxation function in the sol phase away 
from criticality, where a is a non-critical a nd p ossibly 
non-universal exponent. We will return to (|26| ) i n the 
next section. 

From @, @ and @ we deduce *i - e^ for e \ 
with an exponent given by the scaling relation 



£ = k + z 



1 



(27) 



Since A > 0, we have z > k and the scaling relation 
@ is compatible with the inequality (pi). Eq. @ 
was obtained previously in [ p9[ from a model density of 
relaxation times with a sharp upper cut-off. 

According to (Q) the viscosity exponent for the Rouse- 
type model under consideration is given by k rj 0.71, 



VI. TIME-DEPENDENT NORMAL-STRESS 
RESPONSE 

First, let us focus on the normal-stress response to the 
inception of a steady shear flow n(t) = kq Q(t). Here & 
denotes the Heaviside-unit-step function. For each reali- 
sation of the crosslinks Eq. (||) leads to 

Wi(t) := a xx (t) - <r yy (t) = 2k 2 f dt't'x(t') (29) 



in accordance with the principle of frame invariance pd| . 
Eqs. (^|) and (|l^) then imply that for all crosslink con- 
centrations below Ccrit the first normal-stress difference 
increases towards its steady-state value like a stretched 
exponential 



Ni (t) = p K 2 y 1 -2n 2 / dt' t' x{t' 



(30) 



with the same exponent a as the shear-relaxation func- 
tion (p6|). In contrast, for c = c cr ;t we deduce from (piJ) 
the algebraic growth N\(i) ~ t 2 ~ A for long times, a re- 
sult already known on a more phenomenological basis 

Second, we consider a sudden shearing displacement 
K(t) = E5(t), where 5 denotes the Dirac-delta function. 
From (||) we infer 



»«)-£*«<)&{£_** 



(s) 



E\(t), 



(31) 



which, after averaging over disorder, amounts to the 
Kohlrausch decay ( 26 ) in the long-time limit for systems 
below the critical point, respectively to the algebraic de- 



VII. OUTLOOK 



cay t for c = c cr it- Upon comparing (31) to the cor- 



responding result cr xy (t) = E x(t) for shear stress, the 
Lodge-Meissner rule Ni(t)/a xy (t) = E, see e.g. §3.4.e in 
PJ , holds for each crosslink realisation in this Rouse-type 
model. 

Third, we consider the double-step strain flow n(t) — 
E 6(t) — E S(t — ti) with t > ti > 0. In this case one can 
verify in an analogous manner the corresponding relation 
Ni{t)/a xy {t) = —E, which is known |l9| to be valid for 
class I simple fluids. 



We hope to stimulate detailed experimental investiga- 
tions on the crosslink dependence of normal stresses close 
to the gelation transition. If such experimental results 
were at hand, one could judge the effects of the simplifi- 
cations which underly the above Rouse-type model, such 
as the neglect of the excluded-volume interaction and of 
the hydrodynamic interaction. 
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